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We present measurements of the decays B + -> X(3872)K + and B° X(3872)K° with 
X(3872) — ► J/ip-K + -K~. The data sample used, collected with the BABAR detector at the PEP- 
II e + e~ asymmetric-energy storage ring, corresponds to 455 x 10 BB pairs. Branching fraction 
measurements of B{B+ -> X(3872)K+) x B(X(3872) -> J/tpn+Tr') = (8.4 ± 1.5 ± 0.7) x 10~ 6 and 
B(B° X(3872)K°) x S(X(3872) -> J/V-tt+tO = (3.5 ± 1.9 ± 0.4) x 10~ 6 are obtained. We set 
an upper limit on the natural width of the X(3872) of F < 3.3MeV/c 2 at the 90% confidence level. 
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The X(3872) was discovered in 2003 by the Belle Col- 
laboration 1] which reported the observation of a narrow 
resonance in exclusive B^ decays to (</A/> tt + 7T ~). 
The new state was then confirmed by CDF [2j, DO [3j, 
and BABAR Hi]. 

There has been great interest in this narrow state, 
with numerous theoretical interpretations having been 
proposed, including a D°D*° molecule, a diquark- 
antidiquark, a tetraquark state, a hybrid charmonium or 
a classical charmonium state. The diquark-antidiquark 
model pi predicts the A (3872) states to be produced at 
equal rates in B° and B + decays with a mass difference 
of (8 ± 3)MeV/c 2 . The S-wave molecule model pre- 
dicts the neutral B branching fraction to be much smaller 
than the charged B one. 

Studies of angular distributions by CDF [8| favor the 
quantum number assignment J PC = 1 ++ or 2 '". 

The X(3872) has also been observed in the A (3872) -> 
J/fy>7 decay mode by BABAR [9], indicating that it must 
have positive C-parity. Therefore, the 7r + 7r~ pair in the 
A(3872) — > J/ipTt + Tr~ must have a negative C-parity and 
an odd orbital angular momentum L, to satisfy C(7r + 
tt-)=-1=(-1) l + s , with 5 = 0. The tt+ tt~ invariant 
mass distribution has been studied by CDF and 
found to be consistent with a p° meson, where the J/tp 
and the p° are in a relative 5-wave. 

Both BABAR [11] and Belle [I2| have observed the 
A (3872) -> D°D*° decay. These searches were moti- 
vated by the fact that the X(3872) was barely above the 
D°D*° mass threshold. The mass measurement results 
are very consistent between the two experiments but are 
about 3 MeV/c 2 (representing ~ 4 standard deviations) 
above the mass measured in the J/%l)it + ir~ decay mode. 

We report herein an analysis of B + — » A(3872)AT + and 
B° -> A(3872)A°, with A(3872) -► J/^tt+tt- , where 
charge conjugation is implied throughout. We present 
updated branching fractions for the two channels and 
extract the mass and width of the X(3872) state. The 
data sample used for this analysis, collected by the BABAR 
detector at the PEP-II asymmetric-energy e + e~ storage 
ring operated at the Stanford Linear Accelerator Center, 
corresponds to a total integrated luminosity of 413 fb _1 , 
recorded at the T(4S) resonance. 

The BABAR detector is described in detail in Ref. [3] • 
Charged particle momenta are measured with a 5-layer 
double-sided silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) inside a 1.5-T superconducting 
solenoidal magnetic field. A calorimeter consisting of 
6580 CsI(Tl) crystals (EMC) is used to measure elec- 
tromagnetic energy. A ring-imaging Cherenkov detector 
(DIRC) is used to identify charged hadrons, aided by the 
dE/dx measurement in SVT and DCH. Muons are identi- 
fied by the instrumented magnetic flux return (IFR). Par- 
ticle attributes are reconstructed in the laboratory frame 
and then boosted to the e + e~ center-of-mass (CM) frame 
using the asymmetric beam energy information. We use 



a GEANT 4 [3] Monte Carlo simulation (MC) to esti- 
mate the signal efficiencies, employing a sample in which 
one of the generated B mesons decays to the signal mode 
and the other to a representative sample of B decays. 

Charged particles are required to have transverse mo- 
menta greater than 100 MeV/c in the laboratory frame. 
The distance of closest approach of charged tracks must 
be within ±10 cm of the z coordinate (along the beam 
axis) of the primary vertex and within a circle of radius 
1.5 cm in the x—y plane. Kaons, electrons and muons are 
separated from pions based on information from the IFR 
and DIRC, energy loss in the SVT and DCH (dE/dx), 
and the ratio of the associated EMC energy deposition 
{E cs d) to its momentum (E ca \/p). 

The B+ -> J/ip%+%-K+ and B° -> J/tfnr+ir- K° 
decays are reconstructed as follows. Electrons and 
bremsstrahlung photons satisfying 2.95 < m(e + e _ (7)) < 
3.14 GeV/c 2 are used to form J/ip — > e + e~ candi- 
dates. A pair of muons within the mass interval 3.06 < 
m(p + p~) < 3.14 GeV/c 2 is required for a J/i\> — > 
candidate. A mass constraint to the nominal J/ip 
mass [l5| is imposed in the fit of the lepton pairs. We 
reconstruct K® — > ir + ir~ candidates from pairs of oppo- 
sitely charged tracks forming a vertex with a \ 2 proba- 
bility greater than 0.1%, a flight-length (I) significance 
l/cr(l) > 16 (where a{l) is the measurement error) and 
an invariant mass within 15 MeV/c 2 of the nominal 
mass jl5j. We form A (3872) candidates by combining 
J/tp candidates with two oppositely charged pion candi- 
dates, all fitted to a common vertex. Finally, we form 
B + (B°) candidates by combining A (3872) candidates 
with K + {K°) candidates. To suppress continuum back- 
ground, we select only events with a ratio of the second 
to the zeroth Fox- Wolfram moment [l|l] less than 0.5. 

We use two kinematic variables to identify signal 
events coming from B decays: the difference between 
the energy of the B candidate and the beam energy, 
AE = Eg — \fs/2, and the energy-substituted mass 
toes = \J {a/2 + p, • p B ) 2 /E? - p|. Here (E t ,pi) is 
the four-vector (in the laboratory frame) and ^fs is the 
center-of-mass (CM) energy of the e + e~ system, Eg is 
the energy of the B candidate in the CM system and ps 
the momentum in the laboratory frame. The signature 
of signal events is AE ~ 0, and toes_~ m B, where uib 
is the nominal mass of the B meson [l 51 ] . 

If there are multiple candidates in a single event (about 
9% of the events), we select the candidate with the small- 
est value of |A_E|. We optimize the signal selection crite- 
ria by maximizing the ratio rif c / (a/2+^/n™ c ) [13], where 
a = 3 represents the desired significance of signal-to- 
background ratio in number of sigmas and n™° (w™ c ) are 
the number of reconstructed Monte Carlo signal (back- 
ground) events. The optimization is performed by vary- 
ing the selected ranges of AE, |toes — ms], A(3872) 
and K° candidate masses, the K° flight significance and 
the particle identification (PID) selection requirements 
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for the leptons, pions and charged kaons. The criteria 
\AE\ < 20 MeV and \m ES -m B \ < 6 MeV/c 2 , which rep- 
resent about three standard deviations of the resolution 
of the quantities, were found to be optimal for selecting 
signal events. 

We extract the number of signal events with an ex- 
tended, unbinned maximum-likelihood fit to the nix dis- 
tribution, where mx is the J/?pir + Tr~ invariant mass. 

The probability density function (PDF), normalized 
to the total number of events, is V(mx) — J2t n tPt{'mx) 
where n± is the number of events of category t and Vt 
is the associated PDF. We consider only two different 
event categories: signal and combinatorial background 
(which arises mainly from B decays). The signal PDF 
is modelled by a Lorentzian function that describes both 
the natural width and the experimental resolution. We 
model combinatorial background events by a linear func- 
tion in mx- For the neutral mode fit, the width of the 
Lorentzian has been fixed to the value obtained from the 
charged mode fit. 

The fit is performed in the region 3.8 < nix < 
4.0 GeV/c 2 , after applying the optimized selection cri- 
teria on all other variables. The signal region projec- 
tions of the one-dimensional fit to the data are shown 
in Fig. [T] for the B + (top) and B° (bottom) modes. 
We obtain 93.4 ± 17.2 signal events for the B + mode 
(n+) and 9.4 ± 5.2 signal events for the B° mode (n°). 
We interpret the observed events in either mode as the 
X(3872). Results are summarized in Table Q] We fit the 
J/ipir + Tr~ system invariant mass in the toes side band 
region (toes < 5.27 GeV/c 2 ) and observe no signal. 

The efficiency is determined from MC samples with a 
X(3872) signal of zero natural width at 3.872 GeV/c 2 . 
The decay model consists of the sequential isotropic de- 
cays B X{3872)K, X(3872) J/ipp°, and p° -> 
7r + 7r~. This yields a more accurate description of the 
observed ir + ir~ invariant mass distribution [J], compared 
to a three-body decay. Efficiencies are corrected for the 
small differences in K ° reconstruction efficiencies that are 
found by comparing data and MC control samples. The 
final reconstruction efficiencies are (20.60 ± 0.10)% for 
the B+ -> X(3872)K+ mode and (14.50 ± 0.09)% for the 
B° — » X (3872) Jf" mode, where the errors are dominated 
by the size of the signal MC samples. 

The fit is validated with MC experiments, where we 
embed samples of the number of expected signal events 
into MC background samples. On average, the number 
of signal events found is in good agreement with the sig- 
nal sample size. The fit is further validated with a set 
of parameterized MC experiments, based on the signal 
PDF parameters, which return the number of input sig- 
nal events with no significant bias. 

The systematic errors on the branching fraction are 
summarized in Table HU They include uncertainties in 
the number of BB events, secondary branching frac- 
tions [H^, efficiency calculations due to limited MC 
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FIG. 1: Fits to the m(J/ipir + iT~) distributions of (top) 
B+ -> X(3S72)K+ and (bottom) B° -» X(3872)K% can- 
didates drawn from the 413 fb^ 1 sample. The dashed line 
represents the combinatorial background PDF and the solid 
line the sum of background plus the signal PDF. 



statistics, the MC decay model of the X(3872), PID, 
charged particle tracking, Kg reconstruction, background 
modelling (BM) and those arising from fixing the width 
in the B° -> X(3872)K° mode. The production ratio 
of B° and B + mesons in T(45) decays is taken to be 
1.031 ±0.033 0. 

The total fractional errors, 8.8% and 11.7% for the B+ 
and B modes, respectively, are obtained by adding the 
uncertainties in Table HT1 in quadrature. 

The significance is estimated as y— 2 ln(£o/'Cmax) 
where £ max and Co are likelihoods returned by the nom- 
inal fit and by the fit with the signal yield fixed at zero. 
The estimated statistical significance of each signal is 
8.6cr and 2.3a, for the B + and B° modes respectively. 

Using rig and n~ , the efficiencies, the secondary 
branching fractions and the number of BB events, we 
obtain the branching fractions B(B° -> X(3872)K°) x 
B(X -> J/xPtt+tt-) = (3.5 ± 1.9 ± 0.4) x 10~ 6 and 

8.4 ± 



£(^ 



X(3872)if+) x B(X -> J/tPtt+tt-) 



1.5±0.7)xl0 6 . We also calculate a 90% confidence level 
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TABLE I: Fit results for both B+ -> X(3872)iT f and B° -> 
X(3872)A's modes. In the fit to the B° mode, the width is 
fixed to the value obtained from the B + mode. Errors are 
statistical only. The mass measurements are subsequently 
corrected. 





B + -> X(3872)K + 


B u -> X(3872)K% 


Parameters 






mx.flt 


3870.86 ± 0.60 


3868.13 ±1.53 


mx Lorentz Tat 


5.43 ± 1.52 


Fixed to 5.43 


Linear Background 






slope 


-0.30 ±0.04 


-0.28 ± 0.03 


Yields 






signal ns 


93.4 ± 17.2 


9.4 ± 5.2 



TABLE II: Summary of (fractional) systematic uncertainties 
on the branching fraction measurements for both modes. 





B+ -> X(3872)K+ 


B" -> X(3872)Ks 


Tracking 


1.8 


1.4 


K® correction 


n/a 


0.7 


PID 


1.9 


1.4 


MC Model 


1.3 


0.9 


B counting 


1.1 


1.1 


MC statistics 


0.5 


0.6 


Secondary BF 


3.3 


3.3 


BM 


7.5 


4.2 


Fixed width 


n/a 


10.1 


Total Fractional Error 


8.8 


11.7 



(C.L.) upper limit on the neutral branching fraction as 
B(B° -> X(3872)K°) x B(X -> J/ipn+ir-) < 6.0 x lO" 6 
(90%, C.L.). For the ratio of branching fractions, in 
which most of the systematic errors cancel, we obtain 



R(X) 



B(B° -> X(3872)K°) 
B{B+ -> X(3872)K+) 



0.41 ±0.24 ±0.05, 



where the first (second) uncertainty is statisti- 
cal(systematic). Assuming Gaussian errors, we calculate 
the upper limit R{X) < 0.73 at 90% C.L. 

We use the B+ -> ip(2S) K+ and B° -> ip(2S) K° 
decays, with ip(2S) — > J/ipTi + Tr~, in the ip(2S) mass 
region T5j as control modes. We measure the branch- 
ing fractions and obtain the ratio of neutral to charged 
branching fractions R(ip(2S)) = 0.81 ± 0.05 ± 0.01, in 
agreement with the world average of 0.96 ± 0.11. We also 
use these control modes to correct the A" (3872) mass. 
We fit the J/ipTr + Tr~ invariant mass in the ip(2S) and 
AT (3872) region. We correct the AT (3872) mass mea- 
surement, nixfit, b y t he difference between the ip{2S) 
world average mass 15j, m^(2S), and its measured mass, 
"^(2S),fit, which yields m x = mxfit-m^ 2 s),Rt+m^(2S)- 
The result for the B+ mode is (3871.4±0.6±0.1) MeV/c 2 
and (3868.7 ± 1.5 ± 0.4) MeV/c 2 for the B° mode, where 



the first error is the statistical uncertainty on mxfit and 
the second is the uncertainty on m^^2S),fit and rn^(2S)- 
In the neutral mode we have also included an uncer- 
tainty that arises from fixing the width to the value ob- 
tained from the charged mode. The mass difference of 
the A" (3872) states produced in B° and B + decays is 
Am = (2.7 ± 1.6±0.4)MeV/c 2 . 

The natural width Tx of the A" (3872) is obtained using 
the B + mode by subtracting the full width at half max- 
imum of the resolution function measured from Monte 
Carlo r Rcs from the data r fit , T x = r fit - r Res , with 
both the resolution and the natural width of the AT (3872) 
parameterized by a Lorenztian function. We estimate a 
systematic error on the width by comparing the nominal 
value to the value determined when using a two-Gaussian 
resolution function. We determine the natural width of 
the AT(3872) to be (1.1 ± 1.5 ± 0.2) MeV/c 2 , where the 
first uncertainty is statistical and the second systematic. 
From this result we calculate the 90% C.L. upper limit 
on the natural width T x < 3.3 MeV/c 2 . 

In summary, we have performed an updated study of 
the decays B+ -> X(3872)K+ and B° -> X{3872)K° 
with X(3872) — > J/ipir + Ti~ . The branching fraction 
measurements in the B + and B° modes are in good 
agreement with previous results, with comparable or 
better errors. The ratio of the branching fractions is 
R = 0.41 ± 0.24 ± 0.05 and the observed mass differ- 
ence is Am = (2.7 ± 1.6 ± 0.4) MeV/c 2 , consistent with 
either the molecular or diquark-antidiquark model within 
two standard deviations. Finally, we provide an up- 
dated upper limit of the natural width of the AT (3872), 
T x < 3.3 MeV/c 2 (90%, C.L.). 
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